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ABSTRACT 

New observations of CO ( J = 1 — > 0) line emission from M33, using the 25 element BEARS focal 
plane array at the Nobeyama Radio Observatory 45-m telescope, in conjunction with existing maps 
from the BIMA interferometer and the FCRAO 14-m telescope, give the highest resolution (13") and 
most sensitive (<r rms ~ 60 mK) maps to date of the distribution of molecular gas in the central 
5.5 kpc of the galaxy. A new catalog of giant molecular clouds (GMCs) has a completeness limit of 
1.3 x 10 5 Mq. The fraction of molecular gas found in GMCs is a strong function of radius in the galaxy, 
declining from 60% in the center to 20% at galactocentric radius R ga i ~ 4 kpc. Beyond that radius, 
GMCs are nearly absent, although molecular gas exists. Most (90%) of the emission from low mass 
clouds is found within 100 pc projected separation of a GMC. In an annulus 2.1 < R ga i < 4.1 kpc, 
GMC masses follow a power law distribution with index —2.1. Inside that radius, the mass distribution 
is truncated, and clouds more massive than 8 x 10 5 M Q are absent. The cloud mass distribution shows 
no significant difference in the grand design spiral arms versus the interarm region. The CO surface 
brightness ratio for the arm to interarm regions is 1.5, typical of other flocculent galaxies. 
Subject headings: Catalogs — galaxies individual (M33) — ISMxlouds — radio lines:ISM 



1. INTRODUCTION 

The nearby galaxy M33 is an ideal site to study molec- 
ular gas and star formation in the larger context of a 
disk galaxy. Most of our knowledge of star formation 
and molecular clouds is informed by nearby star form- 
ing regions such as Taurus and Orion, but studies of 
star formation across the Milky Way as a whole are 
confused by our perspective from within the Galactic 
disk. In contrast, the M 33 disk face is visible (i ss 52°, 
ICorbelli fe SaluccHl2000l) . so it is easier to study molec- 
ular gas in relationship to other components of the 
galaxy. In addition, M 33 is near enough (D = 840 kpc, 
iFreedman et al.ll2001h that individual molecular clouds 
can be resolved with a large single-dish telescope or mil- 
limeter interferometer. 

Because of the perspective that observations of M33 
offer, the galaxy has been the target of many large- scale 
studies of molecular gas. IWilson fe Scovilig (|1989h ob- 
served the inner region of the galaxy in CO ( J = 1 — ► 0) 
with the NRAO 12-m and followed up with high reso- 
luti on observations using th e Owens Valley interferome- 
ter (jWilson fc Scov ille 1990). The high resolution obser- 
vations were used to study the properties of individual 
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molecular clouds and compare to clouds in the Milky 
Way. A further study (W ilson et al.l Il997h took advan- 
tage of the disk visibility to study variations in the prop- 
erties of molecular gas across a range of galactocentric 
radii. Following improvements in millimeter-wave in- 
strumentation, large-scale surveys of the ga l axy b ecame 
possible, such as that of lEngargiol a et all ([2003, here- 
after EPRB) who used the BIMA millimeter interferom- 
eter to identify al l GMC s across the star forming disk. 
iRosolowskv et al.l (|2003l ) studied a limited area of the 
disk, comprising about 1/3 of the clouds in the BIMA 
survey, at higher resolution. The study determined in- 
dividual cloud pro perties including v i rial m asses, con- 
firmed the work of IWilson fc Scovilld ([1990), and con- 
strained GMC formation mechanisms by comparing the 
predicted and observed values for GMC angular momen- 
tum. Finally, iHever et alJ (|2004L hereafter HCSY) used 
the FCRAO 14-m to complete a single-dish survey of the 
entire galaxy. Comparison with IRAS data determined 
the local star formation rate as a function of the surface 
density of molecular gas. 

The existing observations, while providing a wealth of 
data, also raise additional questions that merit another 
look at the molecular gas in M33. In particular, com- 
parison of the interferometer-only study of EPRB with 
the single-dish survey of HCSY shows that the interfer- 
ometer recovered only 20% of the CO flux in the galaxy! 
Because the typical size scales of molecular clouds in M33 
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should roughly match the synthesized beam of the inter- 
ferometer, such a small fraction of flux recovery is unex- 
pected. What is the nature of the remaining molecular 
gas? One clue comes from the steep mass distribution of 
GMCs in M33. Unlike the Milky Way EPRB found that 
the mass distribution was "bottom heavy" i.e. most of 
the molecular mass was found in clouds near their sen- 
sitivity limit. Thus, it appears that most of the missing 
flux is comprised of low-mass molecular clouds that the 
interferometer could not detect. This prompts several 
questions about the nature of such clouds: Where are 
these low mass molecular clouds located relative to the 
GMCs? Does the fraction of material found in GMCs 
vary over the galaxy? These questions are not easily an- 
swered by the existing observations, because the diffuse, 
extended emission of the lower-mass clouds is not easily 
detected by a radio interferometer, and the existing sin- 
gle dish observations do not have the angular resolution 
to define individual clouds. 

To study the distribution of low mass clouds in more 
detail, we mapped the inner 2 kpc of M33 using the 
BEARS receiver on the 45-m telescope at the Nobeyama 
Radio Observatory (NRO). The NRO 45-m data give 
better sensitivity than previous observations and allow 
detection of a greater fraction of the CO in the central 
region. We present two new maps: a BIMA+FCRAO 
map of the entire galaxy and a BIMA+FCRAO+NRO 
map of the inner 2 kpc. The maps are complementary 
— the first surveys the entire galaxy at 13" (50 pc) reso- 
lution, sufficient to identify individual GMCs, while the 
second has 50% better point source sensitivity and a fac- 
tor of 3.5 better column density sensitivity, though hav- 
ing poorer resolution (20" = 75 pc) and limited coverage. 
Section [2] of this paper describes the new observations at 
NRO, while SJ3] explains the data processing and produc- 
tion of the maps. Section [4] presents the new GMC cata- 
log. The new catalog demonstrates significant variations 
in the mass distribution across the face of the galaxy, as 
described in Sj5] Section [6] summarizes the results. 

2. NOBEYAMA RADIO OBSERVATORY OBSERVATIONS 

The BEARS array on the 45 m telescope at Nobeyama 
Radio Observatory consists of 25 SIS mixer receivers 
arrang ed in a 5 x 5 array with 41"1 between the ele- 
ments (|Sunada et all 120001) . With the 14" beam size at 
115 GHz, the elements are separated by 3 resolution ele- 
ments on the sky. The basic unit of the M33 observations 
consisted of nine separate pointings arranged in a 3 x 3 
grid, each offset from the previous by 14". This strat- 
egy observes a (3.5') 2 region at half-Nyquist sampling. 
To acquire fully sampled data, this observation must be 
repeated three additional times with a 7" offset in each 
direction. Thus, a fully-sampled map consists of 900 in- 
dividual pointings separated by 7"; a half-sampled map 
contains 225 pointings separated by 14". The BEARS 
receiver requires using the 25-element autocorrelator ar- 
ray as a back end. It has a total bandwidth of 512 MHz 
split into 1024 channels, resulting in a nominal velocity 
resolution of 1.3 km s _1 in the 12 CO ( J = 1 — > 0) line. 

The observation strategy divided M33 into (3.5') 2 re- 
gions defined by the footprint of the receiver array on the 
sky. Figure [T] shows the fields observed, and Table[T]gives 
their positions. For each field observed with half-Nyquist 
sampling, we observed the nine individual position offsets 
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Fig. 1. — M33 field locations observed with the NRO 45-m. 
The boundaries of the observed fields established by the 3.5' size 
of the BEAR S recei ver are indicated on an Ha map of the galaxy 
IMassev et~a!l 12006fl . Of the 14 fields shown, the six fields with 
grey borders were observed with full sampling and the remaining 
8 fields were sampled half as densely. A dashed ellipse is drawn at 
R ga l = 2 kpc, illustrating that the survey covers the galaxy out to 
this distance. Properties of the fields are given in Table [T] 

in 20 second scans, then position switched to an absolute 
reference position. The reference positions were chosen 
to be H II regions at large galactocentric radius (R ga i)- 
These regions may contain CO emission, but they were 
always chosen so that the velocity of any emission in the 
reference beam would be much different from the line ve- 
locity in the observation position. In good observing con- 
ditions, two scans could be obtained between reference 
positions. In marginal weather, the number of scans had 
to be decreased to one, and the reference position was 
an offset of 30' in azimuth instead of an absolute posi- 
tion. One pass through the nine positions of a field took 
five to six minutes depending on the settle time for the 
telescope, and typical observations utilized 12 cycles ob- 
taining 180 seconds of integration time at each position. 
This completed a half-Nyquist observation of a field in 90 
minutes. Before and between each field, we updated the 
pointing solution for the telescope using the S40 receiver 
to observe SiO (J = 1 — > 0, v = 0,1) maser emission 
from IRC+30021. In low wind conditions, the pointing 
solution drifted by < 5" over the hour required to ob- 
serve a single field. When the wind velocity was steadily 
above 5 m s , the pointing drift could exceed 10" over 
an hour. 

We observed from 2005 February 9 to February 17 as 
long as M33 was at suitable elevation, excluding time 
lost to maintenance and bad weather. In that time, we 
completed observations of 14 of the (3.5') 2 ; six fields are 
Nyquist sampled and the remaining eight fields are half- 
Nyquist sampled (see Figure [Hand Table Q] for details). 

Since each element of the BEARS receiver consists of 
a double side band SIS receiver, we calibrated the ar- 
ray with observations using a single-sideband receiver to 
obtain an accurate intensity scale. At the latitude of 
the Nobeyama Radio Observatory, M33 transits above 
the elevation limit of the 45-m telescope (80°) render- 
ing it unobservable for 90 minutes in the middle of the 
observations. During that time we calibrated the ef- 
ficiency of each pixel of the BEARS receiver by com- 



GMCs and CO in M33 



3 



TABLE 1 
NRO/BEARS Field Parameters 



Field 



Central Position 3 



Scans b Frac. Retained 



1 

2 
3 
1 
5 
6 
7 
8 
9 
10 
11 
12 
13 
11 



01:33:34.4 +30:34:51.2 
01:33:34.4 +30:38:23.6 
01:33:34.4 +30:41:42.2 
01:33:34.4 +30:45:07.7 



01:33:50.9 
01:33:50.9 
01:33:50.9 
01:33:50.9 
01:34:06.8 
01:34:06.8 
01:34:06.8 
01:34:06.8 



-30:34:51.2 
-30:38:23.6 
-30:41:42.2 
-30:45:14.5 
-30:34:51.2 
-30:38:23.6 
-30:41:42.2 
-30:45:07.7 



01:34:06.8 +30:48:40.0 
01:34:22.8 +30:48:40.0 



108 
108 
108 
108 
108 
382 
324 
323 
112 
436 
108 
108 
320 
163 



0.94 
0.99 
0.97 
0.84 
0.96 
0.71 
0.84 
0.87 
0.95 
0.73 
0.96 
0.97 
0.71 
0.46 



by such signals can be readily identified because the 
rms residual from the b-spline fit is dramatically larger 
than for the remainder of the observations. The small 
=number of spectra for which the rms residual is 4 
Samplmg ti mcs the median were rejected. Table [1] reports the 
Half-Ny4jraction of scans that remain after wind speed and 
Half_N yfeaseline rejection. The acceptable spectra were Hanning 
y^Jj^ySmoothed and down-sampled to a final channel width of 
Half-Ny^.-6 km s _1 before being gridded into a data cube using 
Nyquisfa, Gaussian smoothing kernel with a FWHM of 20". 

Nyquist 

Nyquist3.2. Merging the Interferometer and Single Dish Maps 
Half-Nyq. 

Nyquist The first step in the analysis was to merge the FCRAO 
Half-Ny^HCSY) and BIMA (EPRB) data sets. As noted in §1, 
Half-Nyi ne interferometer map recovers only ~ 20% of the flux 
Nyquisf oun d m the single-dish map, though the interferometer 
map detects nearly all the GMCs in the galaxy. The 



a J2000 b Number of 20 second scans completed by the BEARS array, 
tion of scans retained after automated scan rejection. 

paring observations of NGC 7538 (a 2 ooo = 23 h ll m 45?5 
Aooo = +61°28'09") made with each element of the ar- 
ray to those made of the same position with the single- 
sideband S100 receiver. The intensity scaling factor for 
each BEARS pixel relative to the S100 receiver was the 
ratio of the integrated intensities in the respective obser- 
vations. The BEARS spectra were multiplied by these 
scaling factors averaged over the course of the observing 
run. The scaling factors ranged from 1.7 to 3.6 with a 
median value of 2.4. Daily variation of the scaling factors 
was < 20% for all elements; this variation dominates the 
uncertainty in the overall flux calibration. This scaling 
places BEARS spectra onto the S100 intensity scale, and 
the spectra must be scaled up by a further factor of 2.56 
(= 1/rjmb) to correct for the main beam efficiency. With 
all corrections, the noise level of the final map is 85 mK 
on the T\ scale. 

3. DATA REDUCTION 

3.1. Nobeyama Radio Observatory Data 

Inspection revealed several pathologies in the 
NRO/BEARS data. These are: (1) pointing errors 
induced by high wind speeds, (2) baseline variations, 
and (3) interloper signals in the base band. Because of 
the data volume — 62300 raw spectra — automated 
processing was needed to correct these problems. Scans 
taken when wind velocities exceeded 7.5 m s _1 were 
rejected in their entirety. Spectral baselines were estab- 
lished for every spectrum by a fourth-order b-spline fit 
with break points established every 100 km s^ 1 , ignoring 
regions within 32 channels of the beginning and end of 
the spectrum, within 20 channels of a bad correlator 
section near channel 300, and within 30 km s _1 of the 
H I line-of-sight velocity at observed and reference po- 
sitions in M33. These exclusions prevented the routine 
from fitting a baseline to actual CO emission in M33 1 . 
The baseline fit accounts for low-order baseline ripples 
of arbitrary shape but cannot account for high-order 
"standing wave" patterns in the resulting data caused 
by interloper signals in the base band. Spectra affected 

The H I velocities were determined as a function of posi- 
tion based on the orien tation parameters and rotation curve of 
ICorbe lli & Saluccl pOOOl l: inclination 52°, major axis position an- 
gle 22°. 



Frai-4-m FCRAO dish recovers all flux at the small spa- 
tial frequencies not sampled by the BIMA observations, 
which have minimum baselines < 7 m. Therefore, the 
two maps can be combined to produce a single high- 
resolution, fully-sampled map. The power in the over- 
lapping range of spatial frequencies shows that the flux 
scales of the single-dish and interferometer maps are con- 
sistent within their uncertainties, and no relative flux 
scaling is required. 

To assess systematic uncertainties in the resulting 
map, the data sets were combined using two tech- 
niques: (1 ) image-domain li n ear c ombination as de- 
scri bed by Stanimirovi c et ail (|1999f ) and implemented 
by I Heifer et al.l (l2003f) and (2 1 Fou rier-domain com- 
bination (iBaiaia fc van Albadal 11979( 1 implemented in 
MIRIAD ijSault et all Il995fh The two methods gave 
results indistinguishable within the uncertainties. Fig- 
ure [21 shows the result produced by the linear method, 
which is used throughout this paper. This map will 
be called the "the merged map." Its final resolution is 
13" x 2.03 km s _1 (a projected beam size of 53 pc), and 
its median noise level is 240 mK. The boundaries of the 
FCRAO and BIMA maps do not exactly coincide, but 
the region of overlap (shown in Figure [2|) covers the ma- 
jority of the galaxy to a radius of R ga i ~ 5.5 kpc. 

3.3. Combining All Data 

The NRO 45-m map shows no systematic differences 
from the merged map. Figure [3] compares the two maps 
in the region of overlap after convolving the merged map 
to match the 20" resolution of the NRO map. The noise 
in the NRO map shows greater variation with position 
than in the merged map, and therefore many positive 
and negative outliers are seen in the pixel comparison 
plot (Figure [3l left). Despite that, the concentration of 
points near the locus of equality shows agreement in the 
relative calibration between the two maps. The spectral 
comparison of the emission from the central 7' (1.7 kpc) 
of the galaxy (Figure^ right) highlights the agreement in 
the flux calibration across the observing band. The plot 
also shows that the combination with the BIMA data 
does not significantly affect the amount of flux found in 
the merged data set. The agreement between all three 
data sets is quite good except near Vlsr — 105 km s _1 
where emission in the reference position of the NRO 45- 
m data produces a negative, though irrelevant, feature 
in the data. 
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Fig. 2. — Map of the integrated intensity of CO emission across M33 from the merged BIMA+FCRAO map. Intensity at each position 
is integrated over a 10 channel (20.3 km s — 1 ) window centered on line-of-sight projection of the rotation velocity at that position. To limit 
the visual confusion, only a single red contour is drawn at a level of three times the value of the local noise. Since there are over 2 X 10 4 
independent positions in the map, some of the emission above this contour level will be due to statistical fluctuations alone. Section [4] 
discusses how spectral information can be used to distinguish real features from statistical outliers. The green boundary demarcates the 
border of the region observed with both BIMA and FCRAO. Only in this overlapping region does the map contain information at all spatial 
frequencies. The center of the galaxy at 02000 = l h 33 m 50?8, 52000 = +30° 39' 36'.'7 is marked with a cross. The boundaries of the spiral 
arms adopted in i|5.3l are shown in blue. 




Fig. 3. — Comparisons of the merged BIMA+FCRAO map with the NRO observations, (left) This panel compares the individual 
data values over the shared region between the two maps on an element-by-element basis showing reasonably good agreement within the 
uncertainties of the map. (right) This panel compares the average spectrum that would be derived for each data cube using a Gaussian 
beam with FWHM of 7' (1.7 kpc) centered on the central position of the galaxy. The three spectra are for the merged BIMA+FCRAO 
data cube (blue), the FCRAO data set alone (black), and the NRO data (red). 
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Fig. 4. — Integrated intensity of CO emission in the central re- 
gion of M33 from the combined NRO+BIMA+FCRAO data cube. 
At each position, the data have been integrated over a 20 channel 
(52 km s~ x ) wide window centered on the line-of-sight projection 
of the rotatio n velocity at th at position as derived from the rota- 
tion curve of Corbclli (2003), and the result is shown in negative 
grey scale. The noise varies across the map because of variable 
integration time and weather, and contours show the statistical 
significance based on local estimates of the noise level. Contour 
levels are —4,-2 (dotted), 2, 4. ..14 times the local rms. The center 
of the galaxy is indicated with a cross and is itself devoid of CO. 

Because the two data sets agree well over their com- 
mon area, a final map of the central region of M33 re- 
quires simply averaging the two data cubes. The aver- 
age was weighted by the inverse square of the noise es- 
timate determined for ea ch pixel in th e data cube where 
the Chauvenet criterion (|Tavlorlll997t ) was applied iter- 
atively to remove high-significance outliers (see EPRB 
for details). The final data cube has a typical rms 
noise temperature of 60 mK and a resolution of 20" x 
2.6 km s" 1 , corresponding to a projected resolution of 
81 pc. At this angular resolution, the noise level of of 
the BIMA+FCRAO map is 100 mK, so combining the 
data yields a 40% improvement in the sensitivity. The 
NRO+FCRAO+BIMA map is shown in Figure H it will 
be referred to here as "the combined map." Because 
of variable integration time and atmospheric conditions, 
the noise varies spatially across the map. Figure O shows 
the la rms noise level as a function of position for the 
final data cube. Table compares the properties of the 
various maps. 

4. REVISITING THE GIANT MOLECULAR CLOUDS OF 

M33 

EPRB presented the first complete catalog of GMCs 
in M33 based on the BIMA survey data alone. The 13" 
(50 pc) resolution of the interferometer data is ideal for 
identifying individual GMCs though insufficient for mea- 
suring cloud sizes. The combination of the FCRAO and 
BIMA maps is even better for identifying GMCs because 
the data are not subject to many of the pathologies as- 
sociated with interferometer-only data such as spatial 
filtering and non-linear flux recovery. In addition, the 
sensitivity of the merged map is marginally increased 
over the BIMA data alone because FCRAO and BIMA 
sample some of the same spatial frequencies. As shown 
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Fig. 5. — Map of the 1<t rms noise in the combined 
NRO+BIMA+FCRAO data cube. The noise varies because of 
differences in integration time across the map and changing at- 
mospheric conditions during the observing. Dark areas indicate 
higher noise as shown by the scale bar on the right. Contours are 
drawn at noise levels of 60, 80, and 100 mK 

in Figure [H the molecular gas is distributed in many 
complexes with low surface brightness filaments connect- 
ing the complexes, which are also connected in velocity 
space. The bright CO complex in the north of the map 
contains the most massive molecular cloud in the galaxy 
(Cloud 1 of EPRB). 

We have generated a new catalog of GMCs in M33 from 
the merged data set. The combined map is less suitable 
for the identification of GMCs since the coarser spatial 
and velocity resolution match GMCs poorly resulting in 
a higher completeness limit in the resulting catalog. The 
merged catalog is restricted to the region of overlap be- 
tween the BIMA and FCRAO maps where all spatial in- 
formation is recovered. The catalog was generated using 
a contouring method, as used by EPRB. The first step 
was to identify all regions in the data set with larger 
than a 2.7oy ms detection in two adjacent, independent 
channels. For each such detection, all pixels with signif- 
icance larger than 2a rms and connected to the 2.7a rms 
"core" became part of the cloud candidate. For each 
candidate, we calculated the probability (P) of drawing 
the most significant spectrum from a random distribu- 
tion of noise (see Appendix A of EPRB). We then multi- 
plied by the number of independent measurements in the 
catalog region (N w 2 x 10 6 ) and considered — ln(A^P) 
as the statistic of merit. The final catalog contains all 
candidates with — ln(NP) > 8.0 and with CO velocity 
centroids within 20 km s -1 of the H I velocity centroid 
(|Deul fc van der Hulstl [l987h at that location. The se- 
lection criteria differ from those of EPRB, who included 
high significance detections across the entire bandpass 
and found four clouds at large velocity separation from 
the atomic gas. Subsequent observations of those appar- 
ent clouds showed that they are not real CO emitters 
but rather could be attributed to malfunctions in the 
BIMA correlator (Rosolowsky, Blitz & Engargiola, un- 
published). We have selected a 20 km s -1 separation 
from the H I velocity since this velocity range contains 
90% of the CO flux (Figure [6]) while containing no spu- 
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TABLE 2 
Map Properties 





FCRAO 


BIMA+FCRAO 


BIMA+FCRAO+NRO 




single dish 


"merged" 


"combined" 


Resolution 


45" X 1.0 km s" 1 


13" X 2.03 km s" 1 


20" X 2.6 km s" 1 


Vlsr Coverage (km s _1 ) 


[-340, 0] 


[-340, 0] 


[-300, -49] 


Maximum R ga i (kpc) 


5.5 


5.5 


2.0 


Mass Sensitivity 3 - (M ) 


7.6 x 10 4 


5 x 10 4 


3.5 x 10 4 


Column Density Sensitivity a (cm — 2 ) 


8.2 x 10 19 


7.5 x 10 20 


2.2 x 10 20 


Central Flux b (K km s" 1 ) 


1.48 ±0.02 


1.47 ±0.02 


1.50 ±0.01 



a For a 4<r rmB detection in two adjacent channels using a CO-to-H2 conversion factor of 2 X 
10 20 cm - 2 (K km s — 1 )~ 1 . b Flux averaged over the central 7' of the galaxy (see also Figure |31 right). 




nn f , ; i 1 

10 20 30 40 

|V C0 -V H ,| (km »-) 

Fig. 6. — Cumulative distribution function (CDF) of CO emis- 
sion found within a given velocity separation of the H I centroid, 
averaged over the entire CO map. Since 90% of the CO emission 
is found within 20 km s" 1 (dotted line), we limit our search for 
GMCs over this range. 

rious sources. 

Table [3] gives the new cloud catalog wherein ob- 
jects are named "M33GMC." Clouds are listed in or- 
der of increasing galactocentric radius. The cloud 
masses are based on a CO-to-H2 conversion factor of 
2 x 10 20 cm~ 2 (K km s -1 ) -1 a nd the previously-quot ed 
distance of 840 kpc to M33 (|Freedman et all 120011) 2 . 
The conversion factor is appropriate for the GMCs in 
M33 and does not appear t o vary significantly ove r 
the inner 4 kpc of the galaxy (Rosolows ky et al.l 12003). 
The method for calculating mass differs from that 
of EPRB; masses here ar e based on the method of 
iRosolowskv fc Lerovl (j2006). This method attempts to 
account for the amount of emission not included above 
the 2oy ms cutoff in the data cube by extrapolating the 
emission profile to the K km s _1 intensity contour. Ac- 
counting for this emission roughly doubles the derived 
cloud masses. Cloud major and minor axes and orien- 
tations are also given in Table [3j however, the sizes are 
not corrected for beam convolution and are thus upper 
limits. Deconvolution of the cloud sizes is unstable for 
large beam sizes (50 pc) a nd low signal-to-noise values 
(|Rosolowskv fc Lerovl[2006l) . 

It is possible that the catalog method adopted in this 
study will blend together individual GMCs into a larger 
complex that, at higher resolution, would be identified as 

2 For other distances, scale cloud mass in proportion to D 2 . 



separate clouds. In their study of this problem, EPRB 
found many instances of cloud group s identified as dis- 
tinct in the high-resolution work of iWilson fc Scoviliel 
(| 19901 ) that were merged in the EPRB catalog. The 
same effects must be present in the current catalog and 
will affect studies of the mass distributions below 32J In 
general, appropriately decomposing a blend moves an in- 
dividual cloud to multiple clouds of lower masses, steep- 
ening the mass distribution and reinforcing any trunca- 
tions in the distribution (provided the mass of the orig- 
inal cloud was near the truncation mass). We proceed 
with this caveat in mind but make no effort to sepa- 
rate the catalog objects into smaller units. Decomposing 
blends is inaccurate at th is physical resolution (50 pc; 
IRosolowskv fc Leroy|[2006l ). 

The new catalog contains 149 clouds to an estimated 
completeness limit of 1.3 x 10 5 Mq. The limit is based 
on linearly scaling the completeness limit derived from 
EPRB's detailed examination of the BIMA-only survey 
to the sensitivity limit in the BIMA+FCRAO data. The 
noise level of the merged map does not vary significantly 
in the region where molecular clouds are found, although 
it does increase in the outer regions of the galaxy. The 
original catalog of EPRB contained 144 clouds (omitting 
the spurious high- velocity outliers) with total mass 1.7 x 
10 7 Mq. Total mass of clouds in the new catalog would 
be 2.0 x 10 7 Mq using the same mass formula as EPRB. 
This increase in the mass estimate stems from inclusion 
of the low surface brightness emission that is not detected 
in the interferometer-only map. The larger increase in 
the catalog GMC mass, to 3.9 x 10 7 Mq, comes from 
accounting for the flux below the 2a rms limit. This total 
mass accounts for about 1/3 of all the mass implied by 
emission in the merged map, 1.2 x 10 8 Mq. (This differs 
from the value reported by HCSY in part because we 
include helium and adopt a different distance and CO-to- 
H2 conversion factor.) The remaining 2/3 of the emission 
is from sources that are of smaller mass than GMCs. 
This represents quite a contrast to the inner Milky Way, 
where the GMCs compri se 80% of the molecular mass 
(jWilliams fc McKedli997h . 

The amount of CO emission from low mass clouds can 
be examined by averaging the data in annuli of galac- 
tocentric radius. Figure [7] compares the radial surface 
density profile for GMCs (M > 10 5 Mq) and for all CO 
emission. Within 1 kpc of M33's center, GMCs comprise 
<60% of the total molecular mass, and the fraction de- 
creases rapidly with radius to <10% near R ga i w 4.0 kpc. 
Beyond that radius, GMCs vanish almost entirely with 
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TABLE 3 
Catalog of GMC Properties 



3GMC 


Position 


Vgmc 


Vhi 




Vgmc — Vhi 


-ln(JVP) 


Mass 


(7 V 


Size b 


Arm 




(a2000,i5200o) 


(km s- 1 ) 


(km s- 1 ) 


(kpc) 


(km s- 1 ) 




(10 5 M ) 


(km s- 1 ) 


(pcxpc) 


Classificatic 


1 


01:33:52.4 +30:39:18 


-169.0 


-169.0 


0.18 


0.0 


62.5 


3.1 ± 0.9 


3.5 + 0.9 


82 x 62 (122°) 
68 x 54 (12°) 


Arm 


2 


01:33:49.5 +30:40:03 


-191.8 


-185.3 


0.18 


-6.6 


42.9 


2.1 ± 0.5 


4.5 ± 1.9 


Arm 


3 


01:33:48.1 +30:39:31 


-168.4 


-169.1 


0.21 


0.7 


15.1 


2.5 ± 2.1 


2.6 ± 2.9 


80 x 62 (105°) 
62 x 53 (85°) 


Arm 


4 


01:33:48.3 +30:38:53 


-158.0 


-159.7 


0.23 


1.7 


32.3 


1.0 + 0.5 


2.3 + 0.8 


Arm 


5 


01:33:52.4 +30:40:32 


-209.9 


-201.4 


0.24 


-8.4 


34.0 


1.7 ± 1.1 


2.7 ± 1.2 


67 x 58 (159°) 
82 x 60 (36°) 


Arm 


6 


01:33:49.6 +30:41:08 


-193.6 


-199.5 


0.46 


5.9 


9.3 


1.6 ± 1.5 


2.4 ± 2.0 


Arm 


7 


01:33:56.8 +30:40:08 


-192.6 


-195.0 


0.46 


2.4 


63.9 


3.5 + 1.0 


3.0 + 0.8 


115 x 63 (105°) 
68 x 58 (103°) 
98 x 50 (35°) 
70 x 53 (34°) 
68 x 60 (29°) 


Arm 


8 


01:33:44.6 +30:38:54 


-163.7 


-161.0 


0.47 


-2.7 


12.3 


1.4 + 1.0 


3.1 + 3.6 


Arm 


9 


01:33:51.2 +30:41:29 


-220.3 


-215.4 


0.50 


-4.9 


14.2 


2.5 + 1.3 


3.6 ± 1.8 


Interarm 


10 


01:33:53.5 +30:41:41 


-224.9 


-219.8 


0.53 


-5.1 


34.8 


1.7 + 0.7 


2.7 ± 1.1 


Intcrarm 


11 


01:33:57.0 +30:41:08 


-213.4 


-211.1 


0.54 


-2.2 


20.0 


1.4 ± 1.2 


2.3 + 0.9 


Interarm 


12 


01:33:55.5 +30:41:37 


-221.8 


-216.4 


0.55 


-5.4 


51.0 


4.8 ± 1.1 


6.5 ± 1.7 


107 x 82 (155°) 
60 x 48 (90°) 
62 x 57 (70°) 


Interarm 


13 


01:33:49.9 +30:37:27 


-149.5 


-149.8 


0.57 


0.1 


12.4 


1.1 ± 1.4 


4.9 + 3.2 


Interarm 


14 


01:33:47.9 +30:41:29 


-203.5 


-200.6 


0.63 


-2.9 


32.1 


0.9 ± 0.5 


1.6 ± 1.4 


Arm 


15 


01:33:52.5 +30:42:15 


-214.0 


-216.2 


0.68 


2.2 


21.7 


2.6 ± 1.6 


3.7 ± 2.6 


91 x 52 (172°) 


Interarm 


16 


01:33:43.1 +30:37:31 


-139.5 


-142.8 


0.69 


3.2 


10.1 


0.4 ± 0.3 


1.8 ± 2.5 


60 x 49 (95°) 
63 x 57 (89°) 
93 x 68 (107°) 
165 x 68 (1°) 
66 x 59 (87°) 


Interarm 


17 


01:33:54.8 +30:37:43 


-167.7 


-164.9 


0.69 


-2.8 


88.9 


2.2 ± 0.4 


3.0 ± 1.0 


Arm 


18 


01:33:59.9 +30:40:45 


-209.8 


-207.3 


0.70 


-2.4 


204.3 


7.7 + 0.9 


4.5 ± 0.6 


Interarm 


19 


01:33:52.9 +30:37:17 


-152.7 


-154.0 


0.71 


1.3 


31.5 


4.7+ 1.4 


4.1 ± 1.3 


Interarm 


20 


01:33:45.3 +30:36:52 


-144.2 


-139.9 


0.73 


-4.3 


30.6 


1.9 ± 1.7 


2.4 ± 1.5 


Interarm 


21 


01:33:59.7 +30:41:32 


-221.5 


-217.3 


0.74 


-4.2 


64.2 


4.5 ± 1.1 


4.9 ± 1.2 


82 x 64 (105°) 
75 x 66 (128°) 
81 x 59 (165°) 
73 x 52 (135°) 
80 x 54 (76°) 
137 x 72 (7°) 
73 x 50 (93°) 
70 x 69 (49°) 
86 x 67 (85°) 


Interarm 


22 


01:33:41.0 +30:39:14 


-166.2 


-165.6 


0.77 


-0.6 


192.9 


6.0 ± 0.7 


4.0 ± 0.9 


Arm 


23 


01:33:57.8 +30:42:21 


-233.7 


-229.9 


0.77 


-3.8 


14.3 


1.7+1.0 


3.2 ± 1.6 


Interarm 


24 


01:33:43.5 +30:41:09 


-195.0 


-195.7 


0.83 


0.7 


37.2 


2.2 ± 1.0 


3.9 ± 1.3 


Arm 


25 


01:33:44.3 +30:41:39 


-198.8 


-198.8 


0.88 


-0.1 


15.9 


1.7 ± 1.4 


4.1 ± 2.2 


Arm 


26 


01:33:55.8 +30:43:02 


-227.3 


-226.6 


0.88 


-0.7 


36.5 


6.3 + 1.6 


3.3 ± 1.2 


Intcrarm 


27 


01:33:39.1 +30:38:00 


-149.4 


-148.8 


0.91 


-0.7 


17.5 


1.3 ± 1.2 


3.0 ± 2.2 


Arm 


28 


01:33:44.9 +30:36:01 


-138.7 


-138.5 


0.93 


-0.2 


61.6 


3.6 + 1.2 


4.0 + 0.9 


Intcrarm 


29 


01:34:02.3 +30:39:12 


-201.1 


-195.4 


0.98 


-5.6 


38.3 


4.2 ± 1.3 


4.1 ± 2.5 


Arm 


30 


01:33:57.2 +30:36:41 


-166.7 


-161.9 


1.08 


-4.9 


36.3 


1.1 ± 0.8 


1.7 + 0.7 


65 x 61 (130°) 
173 x 62 (83°) 
56 x 51 (179°) 
79 x 60 (56°) 


Arm 


31 


01:34:02.8 +30:38:38 


-185.8 


-186.4 


1.09 


0.6 


33.6 


5.9 + 2.0 


6.6 + 2.8 


Arm 


32 


01:33:42.1 +30:35:29 


-128.0 


-129.7 


1.11 


1.5 


12.1 


0.6 + 0.5 


2.4 ± 2.3 


Interarm 


33 


01:33:41.3 +30:41:47 


-205.6 


-196.1 


1.11 


-9.5 


16.6 


2.5 ± 0.9 


3.1 ± 1.2 


Arm 


34 


01:33:36.9 +30:39:30 


-167.7 


-163.6 


1.13 


-4.1 


43.1 


6.4 ± 1.5 


3.7 ± 0.7 


138 x 63 (88°) 


Arm 


35 


01:33:41.6 +30:35:14 


-127.2 


-130.5 


1.17 


3.2 


8.1 


2.0 ± 2.5 


5.9 ± 5.0 


64 x 56 (136°) 
69 x 64 (60°) 


Intcrarm 


36 


01:34:01.3 +30:43:55 


-229.2 


-228.3 


1.20 


-0.9 


44.5 


3.0 + 0.8 


2.9 ± 1.0 


Intcrarm 


37 


01:33:35.7 +30:36:25 


-134.0 


-132.8 


1.25 


-1.3 


54.4 


5.0 ± 1.1 


2.7 ± 0.8 


133 x 72 (112°) 
69 x 59 (22°) 
75 x 54 (73°) 
62 x 54 (46°) 


Interarm 


38 


01:33:34.8 +30:37:02 


-136.0 


-137.7 


1.26 


1.6 


18.9 


1.5 + 0.6 


2.3 + 0.6 


Interarm 


39 


01:34:07.2 +30:41:41 


-217.7 


-215.2 


1.27 


-2.5 


8.9 


2.0 ± 1.2 


2.7 ± 1.4 


Interarm 


40 


01:34:06.8 +30:39:28 


-202.5 


-200.4 


1.32 


-2.1 


38.0 


1.3 + 0.6 


4.1 ± 1.8 


Interarm 


41 


01:33:59.2 +30:36:08 


-153.3 


-150.7 


1.33 


-2.6 


85.3 


2.8 ± 0.6 


3.8 ± 1.1 


101 x 57 (141°) 
79 x 67 (21°) 


Arm 


42 


01:34:01.7 +30:36:43 


-153.7 


-154.7 


1.37 


1.0 


40.2 


3.1 ± 1.1 


2.6 + 0.9 


Arm 


43 


01:33:32.1 +30:36:59 


-144.7 


-141.4 


1.45 


-3.3 


14.5 


0.3 ± 0.3 


1.6 ± 1.1 


49 x 48 (142°) 
75 x 54 (101°) 
82 x 64 (131°) 
75 x 61 (85°) 


Interarm 


44 


01:33:31.9 +30:37:58 


-152.5 


-153.0 


1.46 


0.4 


9.6 


1.4 ± 0.9 


3.3 ± 1.4 


Intcrarm 


45 


01:33:48.4 +30:44:44 


-227.4 


-224.7 


1.46 


-2.7 


8.2 


2.1 ± 1.3 


3.3 ± 1.8 


Arm 


46 


01:34:08.8 +30:39:11 


-194.1 


-193.2 


1.51 


-0.9 


119.5 


3.4 + 0.7 


3.0 + 0.8 


Interarm 


47 


01:33:50.4 +30:33:59 


-134.0 


-129.9 


1.52 


-4.2 


30.5 


4.1 ± 0.9 


3.6 ± 1.1 


111 x 59 (54°) 
83 x 62 (79°) 
82 x 57 (37°) 


Arm 


48 


01:34:06.3 +30:37:41 


-163.5 


-164.7 


1.52 


1.2 


51.8 


2.9 + 0.7 


3.5 ± 1.0 


Arm 


49 


01:33:33.5 +30:41:22 


-185.6 


-180.3 


1.63 


-5.4 


40.9 


2.5 ± 0.7 


2.6 ± 1.1 


Interarm 


50 


01:33:42.8 +30:33:10 


-119.3 


-119.6 


1.65 


0.4 


49.9 


2.7+ 1.6 


2.8 ± 1.1 


82 x 68 (6°) 
77 x 62 (70°) 


Arm 


51 


01:34:11.6 +30:43:15 


-220.0 


-219.7 


1.66 


-0.3 


29.9 


2.3 + 0.8 


2.8 ± 1.2 


Interarm 


52 


01:33:37.4 +30:43:11 


-201.0 


-197.2 


1.68 


-3.8 


17.6 


0.9 ± 0.6 


2.3 ± 1.1 


59 x 53 (164°) 
64 x 58 (87°) 


Arm 


53 


01:34:04.0 +30:35:57 


-146.7 


-146.7 


1.69 


-0.1 


14.4 


0.7 + 0.5 


1.6 ± 1.3 


Arm 


54 


01:34:01.1 +30:46:09 


-238.6 


-237.8 


1.69 


-0.7 


29.0 


1.5 ± 1.2 


3.8 ± 2.3 


76 x 55 (30°) 
64 x 53 (96°) 
79 x 61 (113°) 
86 x 52 (42°) 
66 x 50 (46°) 
74 x 66 (108°) 

73 x 59 (79°) 
79 x 63 (29°) 

74 x 56 (51°) 


Arm 


55 


01:34:08.6 +30:37:40 


-156.7 


-168.8 


1.70 


12.2 


9.7 


0.7+1.0 


2.3 ± 1.9 


Arm 


56 


01:33:59.9 +30:34:49 


-144.5 


-140.5 


1.70 


-3.9 


43.3 


2.4 ± 1.1 


3.5 ± 1.6 


Arm 


57 


01:34:12.9 +30:42:00 


-209.8 


-209.5 


1.70 


-0.2 


24.0 


1.8 ± 1.0 


2.7 ± 1.6 


Interarm 


58 


01:33:46.7 +30:45:29 


-232.2 


-227.6 


1.73 


-4.5 


8.2 


1.2 ± 0.9 


2.3 ± 1.8 


Arm 


59 


01:34:10.8 +30:44:53 


-234.6 


-236.0 


1.78 


1.4 


18.9 


2.7+ 1.5 


3.2 ± 2.5 


Interarm 


60 


01:33:52.4 +30:33:10 


-127.6 


-125.0 


1.79 


-2.5 


14.3 


1.9 ± 1.4 


3.2 ± 1.4 


Arm 


61 


01:33:46.9 +30:32:41 


-117.6 


-114.3 


1.80 


-3.3 


73.4 


3.4 + 1.2 


4.1 ± 1.4 


Arm 


62 


01:34:02.6 +30:46:32 


-245.6 


-242.2 


1.80 


-3.4 


34.7 


1.9 ± 0.6 


3.2 ± 1.2 


Arm 


63 


01:33:52.5 +30:46:28 


-240.2 


-238.1 


1.82 


-2.1 


17.5 


1.4 ± 1.3 


3.0 ± 3.2 


69 x 52 (173°) 
72 x 63 (121°) 
69 x 55 (95°) 


Arm 


64 


01:34:04.3 +30:46:33 


-248.8 


-244.2 


1.84 


-4.6 


11.8 


3.0 + 2.0 


2.5 + 2.8 


Arm 


65 


01:33:29.0 +30:40:23 


-179.6 


-177.0 


1.85 


-2.6 


21.3 


1.3 + 0.9 


4.2 ± 2.6 


Intcrarm 


66 


01:34:13.3 +30:39:06 


-196.2 


-196.5 


1.89 


0.3 


15.5 


1.1 ± 0.9 


2.9 + 2.4 


75 x 48 (141°) 
66 X 50 (65 ) 


Interarm 


67 


01:33:38.5 +30:32:20 


— 123.6 


— 119.4 


1.89 


—4.2 


13.9 


0.9 ± 1.4 


2.8 ± 2.0 


Arm 


68 


01:33:44.2 +30:32:05 


-106.3 


-111.9 


1.93 


5.6 


11.7 


0.5 + 0.7 


1.7 + 2.4 


62 x 53 (179°) 
78 x 61 (59°) 


Arm 


69 


01:34:10.8 +30:37:11 


-168.9 


-169.0 


1.95 


0.1 


32.8 


4.1 ± 1.4 


4.4 + 1.7 


Intcrarm 


70 


01:33:37.9 +30:44:45 


-211.0 


-210.3 


2.00 


-0.7 


34.2 


3.7+1.8 


5.3 ± 1.6 


79 x 59 (128°) 
65 x 58 (70°) 


Interarm 


71 


01:33:37.0 +30:31:58 


-121.5 


-116.2 


2.01 


-5.3 


54.0 


0.9 + 0.6 


1.8 + 0.7 


Arm 


72 


01:33:32.1 +30:32:29 


-136.7 


-126.5 


2.03 


-10.2 


16.9 


0.4 + 0.3 


2.0 + 1.3 


56 x 47 (179°) 


Arm 


73 


01:33:33.4 +30:32:14 


-130.4 


-122.3 


2.03 


-8.1 


25.1 


1.0 + 0.8 


2.4 + 1.5 


71 x 52 (107°) 
66 x 58 (115°) 
62 x 57 (61°) 


Arm 


74 


01:34:12.4 +30:37:16 


-189.0 


-179.9 


2.06 


-9.1 


44.3 


1.8 + 1.3 


4.1 ± 1.9 


Interarm 


75 


01:33:33.2 +30:31:56 


-110.3 


-116.9 


2.10 


6.5 


30.5 


0.8 + 0.4 


1.5 + 0.6 


Arm 


76 


01:34:10.7 +30:36:15 


-159.2 


-164.3 


2.11 


5.0 


128.2 


3.5 + 0.5 


1.7 + 0.7 


90 x 64 (168°) 


Interarm 


77 


01:34:16.6 +30:39:17 


-189.5 


-191.1 


2.13 


1.6 


59.6 


3.7 + 0.9 


2.6 ± 1.3 


85 x 67 (47°) 
126 x 59 (73°) 


Interarm 


78 


01:33:30.3 +30:32:16 


-139.9 


-125.1 


2.14 


-14.8 


17.4 


1.7 + 2.7 


2.4 + 2.3 


Arm 
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Fig. 7. — Radial surface density of molecular gas mass. The top 
panel shows the total mass surface density (upper curve) and the 
surface density of mass in GMCs (M > 10 5 Mq — lower curve). 
The bottom panel shows the ratio between these quantities. Error 
bars show the uncertainties in each 500 pc bin. At nearly all radii, 
GMCs comprise a minority of the total molecular gas mass, and 
the fraction declines from the galactic center. There appears to be 
a sharp cutoff in the distribution of GMCs at R ga i ~ 4.0 kpc, but 
the total molecular gas surface density profile continues beyond 
that radius. 

only a few, low-mass clouds beyond this radius although 
there is plenty of molecular gas. The dearth of GMCs 
beyond 4 kpc is real: the catalog completeness limit 
is only 10% higher in mass at R ga i — 5 kpc than at 
Rgai — 0. Thus there are significantly fewer GMCs per 
unit total molecular mass beyond R ga i = 4.1 kpc than 
within that radius. The radius 4 kpc corresponds to the 
outer ends of th e major spiral arms in M33 identified 
in N I and S I (jHumphrevs fe Sandagel 119801 ). In ad- 
dition, the neon abundance shows a sharp drop at this 
radius (|Willner fc Nelson-Pa tel 2002). However, star for- 
mation as traced by Ha continues out to R ga i = 6.7 kpc 
(|Kennicuttlll989f). and the dis k is gravitationally unstable 
to that radius (|Corb clli 2003). Whatever the explanation 
for the changing conditions at R ga i > 4 kpc, star forma- 
tion in the outer part of M33 must occur in molecular 
clouds with M < 1.3 x 10 5 Mq. It is worth noting that 
this effect confounded the total mass estimate of EPRB 
who extrapolated the GMC mass to the total mass of the 
galaxy based on the fraction of emission in GMCs for the 
center of the galaxy. 

The location of the low mass clouds in the inner re- 
gion of M33 is shown directly in the combined map (Fig- 
ure |S]| with its 50% lower rms noise. The size of ob- 
jects identified in the merged map is indicated with el- 
lipses, representing the size and orientation of the ob- 
jects as they would appear in the combined map after 
accounting for the differing resolutions and sensitivities. 
Emission beyond the boundaries of the identified GMCs 
arises from low mass clouds. This emission forms a fila- 
mentary structure, generally surrounding and linking the 
high mass GMCs. The GMCs themselves appear to de- 



45'0' 



40'0' 



30°35'0' 



01"34 m 30 : 




s 33 m 45 s 
a [J2000] 



Fig. 8. — Map of the integrated intensity of CO emission 
in the central region of M33 in gray scale from the combined 
NRO+BIMA+FCRAO map. The gray scale runs linearly from 
to 8 K km s — 1 and a contour is drawn at 1 K km s — 1 . Positions 
of GMCs are indicated with red (blue) ellipses for clouds above 
(below) the 1.3 X 10 5 Mq completeness limit of the revised catalog 
from the BIMA+FCRAO map. The extent of the ellipses is deter- 
mined by the size of the object in the merged map, extrapolated 
to the zero intensity contour and convolved to the resolution of the 
combined map. The border of the NRO observations is shown as 
a green line. Faint CO emission associated with low mass clouds 
is found around the high mass molecular clouds. 

fine the compact peaks of the lower surface brightness 
structures. The low mass clouds are clearly correlated 
with the high mass GMCs, and we quantify this rela- 
tionship by measuring the amount of flux found within 
a given separation from the edge of a catalog GMC. The 
results are shown in Figure [9] which plots the cumula- 
tive distribution of CO flux as a function of projected 
separation from the edge of the catalog GMCs shown in 
Figure [51 The curve rises from ~ 30% at zero separation 
(the amount of flux in GMCs) to over 90% at 100 pc sep- 
aration. The dotted curve indicates the fraction of the 
area found within the same separation, which would be 
the curve that the flux distribution would follow if there 
were no spatial correlation between low mass and high 
mass clouds. The excess above the dotted curve shows 
the spatial correlation of the emission. There is no evi- 
dence for a galaxy-spanning, diffuse molecular gas com- 
ponen t, confirming the conclusions of iRosolowsky et all 

pol . 

5. VARIATIONS IN THE GMC MASS DISTRIBUTION 
5.1. Radial Variation 

The cloud size distribution varies with galactocentric 
radius in a surprising way. Fig. [7] shows the absence 
of GMCs beyond 4.1 kpc, but the most massive clouds 
(M > 8 x 10 5 Mq) are also absent inside a radius of 
2.1 kpc. Indeed all six clouds with M > 8 x 10 5 Mq 
are at 2.1 < R ga i < 2.5 kpc. Evidently something in the 
inner galaxy eliminates very large clouds, and something 
in the outer galaxy eliminates GMCs nearly altogether. 

To examine the cloud distribution more quantitatively, 
we consider two regions divided at R ga i = 2.1 kpc. This 
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Fig. 9. — Cumulative distribution of CO flux as a function of sep- 
aration from the edges of the GMCs shown in Figure [8] using the 
combined BIMA+FCRAO+NRO map. The solid line represents 
the flux distribution and the dotted line shows the cumulative dis- 
tribution of the survey area. The excess of the flux curve over 
the area curve shows the spatial correlation of the CO flux from 
low-mass clouds with the cataloged GMCs. Since the separation 
is measured from the edge of GMCs rather than the center, the 
curves begin with the 33% of the emission found in GMCs which 
are found in the 10% of the map area. 
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Fig. 10. — Cumulative mass distribution of GMCs in the inner 
(blue, Rg a i < 2.1 kpc) and outer (black, R ga i > 2.1 kpc) regions of 
M33. A power law cumulative mass distribution is shown for the 
outer galaxy (black dashed line), and truncated and non-truncated 
power law distributions are shown for the inner galaxy (blue dashed 
and red solid lines respectively). The vertical dotted line indicates 
the completeness limit for the survey. 

division yields approximately equal GMC mass in each 
region. Figure [TO] shows the mass functions for the GMCs 
in these two regions. The mass distributions are signifi- 
cantly different. In the outer GMC region (the annulus 
2.1 < R ga i < 4.1 kpc), the GMC mass distribution fol- 
lows a power law from high mass down to the complete- 
ness limit of the catalog. In the inner region, however, 
despite the existence of many more GMCs (about three 
times more per unit area), the mass distribution is trun- 
cated at the high end and may also fall slightly below a 
pow er law distribution at the low end. A two-sided KS 
test ([Press et al.lll992j ) shows there is only a 0.8% proba- 
bility that the two samples could be drawn from a single 
distribution while exhibiting such a large difference. 



The cumulative mass distribution functions can be 
characterized quantitatively by truncated power-law 
functions with the form: 



N(M' > M) = N 



M 
Mo 



7+1 



- 1 



(1) 



where Nq, Mo, and 7 are parameters, and the fit is 
restricted to clouds above the completeness limit of 
1.3 x 10 5 Mq. This functional form is adopted for com- 
parison with the GMC population in the Milky Way. The 
parameter 7 is the index of the mass distribution. If 
7 < — 2, most of the mass is found in low mass clouds, 
and if 7 > —2, the high mass clouds dominate as is the 
case for the inner Milky Way. As long as N is signif- 
icantly larger than unity, the mass distribution has a 
truncation near M — 2 1 '' 7+1 ^Mo. In the inner Milky 
Way, the mass distribution is commonly expressed as 
a differential mass distribution. The differential mass 
distribution is well-characterized by a power-law with a 
truncation at large clou d masses (Mo ~ 6 x 10 6 M Q 
Williams & McKcc 1997, and references therein): 



dN 
dM 



oc 



M 
Mo 



M < Mq. 



(2) 



A similar form of the distribution, albeit with different 
slopes, has been r eported in t he outer Milky Wa y (e.g. 
iHever et alJl200l . the LMC dFukui et aUl200lh . M33 
(EPRB), and other systems farther afield. Fitting the 
cumulative distribution instead of the differential one is 
a better way to describe the data because the cumulative 
function is not affected by biases introduced by binning, 
and it can accoun t for uncertainties in the cloud masses 
(lR,osolowskvl l2b05V Table H gives the best-fit coefficients 
for the M33 data, and the resulting distributions are plot- 
ted Figure [TU1 As noted above, there is no significant 
truncation of the mass distribution for the outer galaxy 
clouds, and Tablc|4]reports coefficients for a simple power 
law fit. For comparison, Table 0] also includes the fit to 
the cumulative mass distribution of all the clouds in the 
catalog. Both the fits and Figure [10] indicate that the 
primary difference between the two distributions is the 
truncation in the upper end of the mass distribution of 
the inner galaxy clouds. No such truncation is apparent 
in the outer portion of the galaxy, and including one in 
the fit does not appreciably change the derived values or 
goodness of fit. In the mass range where the inner galaxy 
distribution fits a power law, 1.3 < M/(10 5 Mq) < 5.0, 
the two distributions are indistinguishable. While the 
inner galaxy has a larger fraction of its molecular gas 
in GMCs (Figure EJ, the formation or survival of the 
highest-mass GMCs (>7x 10 5 M ) is suppressed. One 
possibility is that high mass clouds are sheared apart by 
galactic tides in this region; but this possibility is un- 
likely because clouds with masses of 10 7 M w w o uld b e 
stable 3 according to the criteria of lStark fc Blitzl (|1978h . 

The mass truncation and differences in the slope of the 
distribution may be effected by several aspects of the ob- 
servations and the cataloging procedure, but none of the 
likely effects would produce a truncation where none is 

3 We have assumed the surface densities and virial parameters 
would comparable to ot her GMCs in M33 and used average values 
of these quantities from Blitz ct al. (2007). 
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TABLE 4 

Parameters of Mass Distributions in M33 



Objects 


7 ^0 


M>/(10 5 M Q ) 


All Clouds 


-2.0 ±0.2 11 ±2 


14 ±7 


Radial Variation 


Outer galaxy 
Inner galaxy 
Inner galaxy a 


-2.1 ±0.1 

-1.8 ±0.2 21 ±11 
-2.4 ±0.2 


34 ± 14 
7.4 ±0.5 
23 ±5 


North-South Variation 


North 
North 
South 
South a 


-2.0 ±0.2 5.4 ±4.3 
-2.25 ±0.15 
-1.4 ±0.2 41 ±11 
-2.4 ±0.3 


34 ± 12 
34 ±11 
6.4 ±0.5 
20 ± 16 


Arm-Interarm Variation 


Arm 
Interarm 


-1.9 ±0.2 5±4 
-2.2 ±0.2 4±4 


16 ±5 
13 ±4 



a A mass spectrum wit hout a truncation is also fit to the 
data as shown in Figure HOI to illustrate the significance of 
the cutoff found in the southern/inner galaxy clouds. 



present. IWilson fc Scov ille (1990j) have already reported 
a truncation in the mass distribution of GMCs, but it is 
difficult to assess their results as their survey is spatially 
incomplete and selected targets based on dust and CO 
priors. Since the merged data span the entire galaxy, 
spatial biases will not af fect the mass distribution. The 
IWilson fc Scoville] (|1990f ) data had significantly higher 
resolution and several of the clouds considered as sepa- 
rate in their study are not resolved into individual clouds 
here (or in EPRB, see their §3.5 and Figures 13 and 14). 
Blending will act to minimize the effects of truncation 
since separate small clouds will appear as a single cloud 
of higher mass, so the truncations are certainly present 
in the GMC mass distributions. The correction for emis- 
sion found below the 2<r boundary of the clouds will also 
change the mass distribution slightly. Since flux loss pri- 
marily affects clouds near the completeness limit (since 
most of the emission in these clouds is below the identifi- 
cation contour), correcting for the effects of the flux loss 
will move these clouds to higher mass, thereby steepen- 
ing the spectrum. Since this is a correction for an in- 
strumental effect, making this correction likely improves 
estimates of the mass distribution. Finally, regardless of 
the precise nature of the systematic effects on the mea- 
sured mass distributions, this reported measurement is 
differential between two regions cataloged with a uniform 
method. As such, the reported differences must reflect 
some underlying difference in the character of molecular 
gas in the two regions. 

Similar variations in the GMC mass distribution as a 
function of position in a galaxy have not been directly 
quantified in previous work, but the re is some eyidence 
for such variations in the Milky Way. iRosolowskvl (|2005t ) 
argued for different mass distribution functions between 
the inner and outer Galaxy. This variation within the 
Milky Way is uncertain, however, because differences 
in cloud cataloging methods can produce apparent dif- 
ferences in the mass distribution while the underlying 
cloud populations are the same. In addition, large arm- 
interarm contrasts observed in many spiral galaxies sug- 



gest that the mass distribution of GMCs changes az- 
imuthally (see also §5.3p . 

5.2. North-South Variation 

The right-hand panel of Figure [3] implies that the ap- 
proaching and receding sides of the galaxy have a signif- 
icant flux asymmetry. The northern (approaching) por- 
tion of the galaxy has a total flux of 8200 Jy km s~ 
whereas the southern (receding) half of the galaxy has a 
flux of 7800 Jy km s _ \ a difference of 5%. The asym- 
metry is also seen in the H I (where the approaching 
half of the galaxy has 13% more flux than the receding) 
and the Ha (5%) distributions (based on the images of 
iMassev et al.ll2006D . making M33 slightly lopsided. Moti- 
vated by this flux asymmetry, we also consider the varia- 
tion in the GMC mass distributions between the northern 
and southern halves of the galaxy. We define the south- 
ern (receding) portion of the galaxy as the portion of the 
disk with galactocentric azimuth within 90° of the kine- 
matic major axis and the northern region as the comple- 
ment of the southern region. We fit truncated power-laws 
to the two distributions and report the derived param- 
eters in Table |4] Like the variation between the inner 
and the outer galaxies, there is substantial difference be- 
tween the northern and southern mass distribution, par- 
ticularly at high masses. The southern portion of the 
galaxy lacks high mass clouds that are present in the 
northern portion of the galaxy. Both mass distributions 
have some evidence for truncation at the high masses, 
though the evidence is relatively weak for clouds in the 
northern portion of the galaxy. While it would be inter- 
esting to separate the variation in the mass distributions 
into 4 separate regions (inner vs. outer and northern vs. 
southern), we lack sufficient numbers of clouds for reli- 
able fits. Since differences between the inner and outer 
regions of the galaxy are more pronounced, we focus our 
analysis on these difference, though with better data, the 
north-south asymmetry could be explored. 

5.3. Arm-Interarm Variation 

For galaxies with well defined spiral structure, the 
spiral arms contai n nearly all of the molecular mate- 
rial in the galaxy (fVogel et alj 119991 : iRegan et all 120011 ; 
iNieten e~a l. 2006) . The Milky Way appears to b e such a 
galaxy ( e.g. iDame et alj2001l;lstark fe Led200(l . Molec- 
ular gas in the interarm regions of spiral galaxies appears 
to have different properties than that found in the arms 
(I Walsh et al.l[200alTosaki est alll200l iStark fc Ledl2006l) 
though iTosaki et al.l ( 20031 ) suggest that the properties 
of individual GMCs may not vary significantly. While 
our observations lack the resolution to assess whether 
the GMC properties are substantially different in inter- 
arm regions of M33, we can investigate whether there 
is variation in the mass distributions of GMCs in and 
out of spiral arms. Unfortunately, the spiral arms in 
M33 are not nearly as well defined a s they are in the 
Milky Way. IHumphrevs fc Sandagd (| 1980( 1 identified 
ten spiral arms in the flocc ulent structure of M33 w ith 
two arms being dominant. IHumphrevs fc Sandagd at- 
tributed the two primary arms to spiral density waves. 
This conjecture was bolstered by observations in the 
near infrared, which revealed that the spiral structure 
of the old stellar population, and hence the stellar 
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Fig. 11. — Comparison of the cumulative mass distributions for 
arm and interarm GMCs. Truncated power law fits reveal a slightly 
shallower underlying slope for the "arm" distribution, though the 
level of truncation in the distributions appears similar between the 
two populations. 

mass, was dominated by these two "grand design" spi- 
ral arms (|Reean fe Vogell[l99l . The remaining eight 
arms are mainly located in the outer galaxy. They 
were primarily defined by the presence of young stel- 
lar associations but also correlate well with narrow 
filament s of atomic gas (iDeul fc van der Hulstl H987L 
EPRB). IHumphrevs fc Sandagel attributed these eight 
secondary arms to an unknown mechanism other than 
density waves. 

For present purposes, we have defined the locations 
of the M33 spiral arms using a 3.6 /jm image of the 
galaxy taken usi ng the IRAC instr ument of the Spitzer 
Space Telescope (|Gehrz et al.ll2005T) . A small correction 
was made for dust emission in the 3.6 /im image leaving 
emission primarily from the old stellar population. An 
azimuthally-symmctric exponential disk was fit to the 
result. The spiral arms become apparent as regions of 
positive surface brightness after subtracting the model 
disk. Detailed analysis will be presented elsewhere, but 
Figure shows the spiral arm locations with respect to 
the CO emission. 

Figure [TT1 compares the mass distributions of GMCs in 
and out of the grand design spiral arms of M33. Only 
slight differences are seen. The distributions in both arm 
and interarm regions show a truncation at high mass, but 
massive GMCs are found in both environments. The 
truncation levels and masses are roughly equal. Al- 
though a two-sided KS test suggests that the distribu- 
tions have different shapes, the difference is barely sig- 
nificant (Pks = 0.04). Unlike the sharp cutoff in the 
inner galaxy mass distribution, the arm and interarm 
mass distributions lack a specific feature that sets them 
apart. Table [4] gives parameters of truncated power law 
fits. 

In addition to comparing arm-interarm mass distribu- 
tions of GMCs, we examined arm-interarm variations in 
the column density of molecular material. Figure [T2] 
shows the ratio of the average surface density of CO 
found in spiral arms compared to that in interarm regions 
as a function of galactocentric radius. In the inner part 
of the galaxy, the ratio has a roughly constant value of 
1.5, comparable to oth er fiocculent galaxies (NGC 5055: 
2-3. iKuno et~all Il997l : NGC 6946: - 1.4. iWalsh et~all 
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Fig. 12. — The ratio of surface densities of molecular gas between 
the arm and interarm regions of M33. The inner regions of the 
galaxy show a constant value of ~ 1.5 but there is an apparent 
increase at larger galactocentric radius. 

2002). The value is significantly smaller than galax- 
ies with stronger grand-design structure such as M 31, 
where the ratio approaches 20 (Ni eten et al.l I2006T ). or 
M51, where contrasts of a factor of 17 have been ob- 
served (|Garcia-Burillo et al.l 1993). In M33, much of the 
"interarm" molecular material is associated with the sec- 
ondar y spiral arms identified bv IHumphrevs fc Sanda"gel 
(1980). However this is a natural consequence of star 
formation associated with GMCs. Wherever GMCs ex- 
ist, one would expect newly formed stars to exist as well 
and create secondary "arms." The strong correlation be- 
tween GMCs and high mass stars (as traced by H II re- 
gions) was established in EPRB. Regardless of the visual 
appearance, in the inner regions of M33, the CO distribu- 
tion is only weakly coupled to the stellar mass distribu- 
tion, as evidenced by the small arm/interarm variation. 

The increase in the arm-interarm contrast at large 
galactocentric radius appears statistically significant. 
However, at these distances, the spiral density wave be- 
comes poorly defined and may be following enhance- 
ments in the stellar light distribution associated with 
star formation in molecular clouds thereby increasing the 
arm-interarm contrast. Carefully determining the loca- 
tion and amplitude of the spiral pattern present in the 
old stellar population will clarify these results. 

5.4. Understanding Variations in GMC Mass 
Distributions 

The results of §ij4] and [5] demonstrate that the mass 
distribution of GMCs is changing significantly across the 
face of the galaxy. The highest mass GMCs are found 
at intermediate galactic radii (2 kpc < R ga i < 4 kpc) 
and no clouds with masses larger than 8 x 10 5 M Q oc- 
cur outside this region. These changes are circumstan- 
tially related to variations in the galaxy roughly demar- 
cated at these radii (^2 and ~4 kpc). To generally de- 
scribe the galaxy, we divide it into three sections at small 
{Rgai < 2 kpc), intermediate (2 < R ga i < 4 kpc) and 
large (R ga i > 4 kpc) galactic radii. At small galacto- 
centric radii, the spiral structur e is not as well-defined 
as it is for intermediate r adii (Humphreys & Sandagc 
Il980t iRegan fc Vogeilll994h . At large radii, the disk of 
the galaxy develops a warp (|Corbelli fc Schneiderl 119971 : 
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iCorbelli fc Saluccill2000f) and ult imately becom es gravi- 
tationally stable (R ga i > 6 kpc: ICorbelilll2003l ). These 
divisions between these gross features of the galaxy cor- 
respond roughly to the regions where we see changes in 
the molecular cloud mass distribution. 

While this correlation is suggestive, the likely regu- 
lator of the molecular cloud masses is structure of the 
atomic gas from which the GMCs must have formed 
(IRosolowskv et al.l [20031: [Blitz et all 120071) . The GMCs 
are found at the peaks of the atomic gas distribution 
(EPRB), a generic feat ure of the ISM in galaxies dom- 
inated by atomic gas (jBlitz et al.l l2007f ). In M33, the 
atomic gas is d istributed in what appears to be a filamen- 
tary network (jDeul fc van der Hulstl Il987t ). The char- 
acter of this network varies with radius in the galaxy 
consisting of small, relatively faint filaments in the in- 
ner galaxy whereas there are large (kiloparsec scale) fila- 
ments in at intermediate radii associated with the spiral 
arms and then large, less bright filaments are large ra- 
dius. We quantify this description by decomposing the 
bright filamentary network into a set of clouds and ex- 
amining the "cloud" properties. While we adopt a simi- 
lar procedure to the segmentation that identifies GMCs, 
such "clouds" do not necessarily represent distinct phys- 
ical objects. The segmentation is a way to measure the 
characteristic sizes and masses of objects in the atomic 
ISM of M33 across the galaxy. 

To generate a catalog of H I objects, we apply a three- 
tie red brightness cut at T a = 64 , 48, and 32 K to the map 
of iDeul fc van der Hulstl (|1987l ). Connected regions of 
pixels are identified at each brightness cut. New regions 
at a given brightness cut are identified as new clouds. 
Regions containing pixels from higher levels have those 
pixels assigned to the closest predefined region. Here 
"closest" means the region that was identified at a higher 
level to which a pixel is connected by the shortest path 
contained entirely within the new region. The algo rithm 
is similar to CLUMPFIND (j Williams et al.lll99l . but 
the coarse contouring levels prevent every local maximum 
from being identified CIS (IS 3b cloud. The shortest-path 
criterion is necessary for decomposing the H I emission 
since the "clouds" blend into a continuous network at the 
low brightness thresholds. Figure [8] suggest that a similar 
criterion would be necessary for identifying GMCs in a 
CO map with higher sensitivity. We reiterate that the 
decomposition is only a means to characterize the sizes 
of structures in the atomic gas across the galaxy in a 
uniform manner. 

In Figure ITBl we show the integrated intensity map 
of IDeul fc van der Hulstl (|1987l ). restricted to the GMC 
survey area, with the ellipses indicating measurements 
of the major (£ m aj) and minor (£ m in) axes of the clouds 
in blue. For each H I cloud of mass M, we calculate 
the characteristic mass (M c har) of objects that would 
form by gravitatio nal instability alon g the filament using 
the prescription of Elmegreen (1994). A major result of 
Elmegreen's analysis including both spiral potential and 
magnetic fields is that the fragmentation of neutral gas 
into "superclouds" (which would be substructures in the 
"clouds" we catalog) can be evaluated in terms of the 
local quantities that characterize the gas, rather than 
using disk-averaged values. The characteristic mass of 
these superclouds is thought to establish the upper cut- 
off observed in the molecular cloud mass distribution in 
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Fig. 1 3. — Map of the integrated i ntensity of H I 21-cm emis- 
sion from Dcul & van der Hulst ( 1987) with locations of H I clouds 
indicated. The clouds are identified based on the algorithm de- 
scribed in i|5.4l The clouds are simply convenient descriptors of 
the filamentary structure of the atomic gas and do not necessarily 
correspond to physically relevant objects. 

the Milky Way (|Elmegreenlll994 iKim et al.ll2003l ). In 
this picture, disk instabilities form atomic superclouds 
which are (marginally) self-gravitating and GMC com- 
plexes form by turbulent fragmentation and collapse in- 
side the superclouds. The total mass of the parent su- 
percloud would then establish the maximum mass of the 
distribution. 

We apply the local fragmentation formalism to all 
the H I clouds emphasizing the principal result of the 
Elmegreen analysis: the use of local conditions to eval- 
uate the characteristic fragmentation masses. This ap- 
plication will illustrate how variations in the cataloged 
atomic cloud properties would translate into variations 
in the characteristic mass of superclouds and thus into 
the cutoffs in the GMC mass distribution. In terms of 
the properties of the identified atomic gas clouds, the 
charac teristic mass of superclouds is given by Elmegreen 
(fl99l as: 

where \x is the mass per unit length of the clouds 
(M/£ ma j) and c g is the three-dimensional velocity dis- 
persion of the gas, taken to be 18 km s _1 for M33 
(jCorbelli fc Schneiderlll997t ) assuming isotropy. We plot 
the characteristic mass of supercloud formation as a 
function of galactocentric radius in Figure [TU The fig- 
ure shows that the atomic structures that could sup- 
port fragmentation into large mass superclouds are most 
abundant at intermediate values of R ga i- Indeed, the 
factor-of-three variation in the maximum mass is consis- 
tent with the idea that the mass of fragmentation con- 
trols the variations in the maximum mass of molecular 
clouds seen between inner and intermediate radii. How- 
ever, the results do not adequately explain the abrupt 
fall-off in the number of GMCs seen at large R ga i since 
there are apparently massive clouds at this radius which 
fulfill the instability criteria. Other additional factors 
may contribute to the dearth of GMCs at this radius. 
The surface density of the H I disk becomes compa- 
rable to the stellar disk at R ga i ~ 4.5 kpc which will 
result in a significant inc rease in the H I scale height 
(|Blitz fc R osolows kvl2006l) , not accounted for in our sim- 
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Fig. 14. — Characteristic mass of fragme ntat ion as a function of 
Rgal for the H I clouds identified in Figure IT3l The mass is calcu- 
lated using the formula from Elmcgrccn (1994) and the measured 
H I cloud properties. The characteristic mass of fragmentation 
peaks at intermediate R ga i which may explain the abundance of 
high mass clouds in that region. 

plincatio n of lElmegreenl (11994 ) . Th e disk also develops 
a wa rp HCorbelli fc Schneider! H997t iCorbelli fe Saluccl 
2000) at this radius though the disk remains gravitation- 
ally unstable t o structure formation out to larger R ga i 
(|Corbellill2003h . 

The locations and mass distributions of GMCs are 
consistent with being governed by the structure in the 
atomic gas. For R ga i < 2 kpc it appears that mas- 
sive stars establish the structure of the ISM based on 
the pre ponderance of wind-driv en shells in the inner 
galaxy (jDeul fe den Hartodfl990l ). The GMCs in the in- 
ner galaxy appear on the edges of these holes (EPRB) 
suggesting that the clouds are formed in the shells swept 
up in winds from high mass stars. At intermediate and 
large (i.e. R ga i > 4 kpc) galactic r adii, there are sub- 
stant ially fewer wind-driven shells (jDeul fc den Hartod 
119901 EPRB), and the molecular ISM becomes more con- 
centrated along the spiral arms (Figure fl"2"|) . At large 
radii, GMCs disappear entirely, as does significant spi- 
ral s tructure, and the rotation curve flattens apprecia- 
bly (jCorbe lli 2003). The molecular gas that is present 
at large radii may be the result of small scale compres- 
sions of the neutral ISM producing less massive molecular 
clouds. Given that applying simple instability analysis 
to real galactic disks only leads to suggestive correla- 
tions, resolving the variation of cloud mass distributions 
will likely require sophisticated simulations of molecular 
cloud formation across the galactic disk with sufficient 
detail to match to these observations. Further progress 
on the effects of spiral structure and molecular cloud for- 
mation will be possible in conjunction the high resolution 
map VLA+GBT map of H I emission by Thilker & Braun 
(forthcoming) . Progress on simulations and observations 
will facilitate the critical evaluation of the above conjec- 
tures and better establish the association between the 
molecular and atomic phases of the ISM. The variations 
in the mass distributions of molecular clouds reported 



in this paper will represent an important observational 
feature to reproduce in theoretical analyses. 

6. SUMMARY AND CONCLUSIONS 

New maps of CO (J = 1 — » 0) emission in M33, pro- 
duced by combining existing interferometric and single- 
dish maps, including new single-dish observations from 
the NRO 45-m, give an improved catalog of giant molec- 
ular clouds (GMCs) and a more sensitive measurement 
of the non-GMC CO emission. The GMC properties in 
the new catalog are unaffected by the spatial filtering 
and non-linear flux recovery of the existing interferom- 
eter catalog. The relative calibrations of all three data 
sets agree, and the total flux in the merged and com- 
bined maps match the FCRAO map in the regions of 
overlap. The GMC properties in the new catalog are 
corrected for instrumental effects using the method of 
iRosolowskv fc Lerovl (|2006f ) resulting in superior esti- 
mates of the cloud properties. 

Improvements in the estimates of cloud properties as 
well as a high sensitivity map of the inner galaxy yield 
several new results. The fraction of molecular gas found 
in GMCs (M > 1.3 x 10 5 M Q ) declines with galacto- 
centric radius, ranging from 60% at R ga i — to 20% 
at R ga i — 4.0 kpc. Based on the higher sensitivity 
BIMA+FCRAO+NRO map, molecular gas that is not 
associated with GMCs is located in diffuse and filamen- 
tary structures around the GMCs. Nearly 90% of the 
diffuse emission is found within 100 pc projected separa- 
tion of a catalog GMC boundary. Beyond 4.0 kpc, the 
GMC mass fraction cuts off sharply, though molecular 
gas is detectable to the edge of the surveyed region at 
Rgai — 5.5 kpc. Even within the area R ga i < 4 kpc, 
where GMCs are found, the mass distribution in the in- 
ner galaxy (R ga i < 2.1 kpc) is significantly different from 
the mass distribution in the outer galaxy. In the annu- 
lus 2.1 < R ga i < 4.1 kpc, the GMCs have a power-law 
mass distribution of dN/dM oc M~ 21 . Inside 2.1 kpc, 
the GMC mass distribution is truncated at high mass 
(> 4 x 10 5 Mq), but lower mass clouds appear to follow 
the same distribution in both regions. We argue that 
variations in the H I gas structure are responsible for the 
observed changes in the GMC mass distribution. 

The southern (receding) half of the galaxy shows 
a truncated mass distribution when compared to the 
northern half of the galaxy. There is no obvious dif- 
ference between the mass distribution of GMCs associ- 
ated with the grand design spiral arms of the galaxy and 
those GMCs in the interarm region. In molecular sur- 
face brightness, M33 shows a arm-interarm contrast of 
1.5, typical of other fiocculent galaxies. 
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